The Kronig-Penney-Ising picture of the colossal magnetoresistance 
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From general arguments, it is shown that a magnetic 
Kronig-Penney model based on the thermodynamics of an 
Ising model can be used for describing the Colossal Mag- 
netoresistance (CMR) phenomenon. The model considers 
a tunneling-like transmission process of hopping electrons 
through a dynamic lattice characterized by evolving magnetic 
clusters. In this model, correlations between the magnetic 
states are considered to be more relevant than the lattice 
strain effects for obtaining the CMR features. Physical ar- 
guments lead to the theoretical description of the intrinsic 
temperature and field dependences of the CMR observed in 
typical manganite materials. 

PACS numbers: 05.50.+q, 75.70.Pa 



It is often argued that the colossal magnetoresistance 
(CMR) is due to a complicated interplay between elec- 
tronic and lattice degrees of freedom (polaron). In the 
following, we show that the magnetic degrees of freedom 
contain the basic contribution to the field and tempera- 
ture CMR behavior. 

Let us recall the intriguing CMR phenomenon. The re- 
sistivity of perovskite materials has an "enormous" peak 
at some temperature Tm-i considered as related to a 
metal-insulator (M-I) transition. This transition is often 
accompanied by a second-order ferromagnetic transition 
at a temperature T c . Under a magnetic field, the p(T, H) 
appreciably decreases near the temperature where p(T, 0) 
has its maximum. Many theories have been proposed for 
CMR, but the exact temperature and field behaviors are 
not yet fully satisfactory understood. 

The CMR phenomenon in perovskite based materials 
is herebelow explained through critical phenomenon be- 
haviors and analytical laws for describing the main ob- 
servations. The model seems much more simple than 
previously presented models. At this stage, it should not 
be requested that it solves all known puzzles for all com- 
pounds with CMR. Simplifications are made herebelow 
for presenting the approach. In fact, we present an al- 
ternative view with respect to the mpst popular theories 
based on the magnetic polaron. idead including the Dou- 
ble Exchange Scenario (DES)I3. Various features seem 
hard indeed to be put into this polaronic-DES frame- 
work, e.g. at high temperature the small polaron pic- 
ture witlua few meV activation energy due to lattice 
distortionu does not directly match to the large polaron 
picture of the metallic state at low temperature, where a 
band picture should emerge and hold below the magnetic 
transitionQ. In fact, the question can be raised whether 



there is a change in the mobility of carriers or in their 
number (or both) or whether there is even a change in 
the band structure necessarily implying some conductiv- 
ity transition. For instance, it can be accepted that the 
mean free path is small at high temperature, but at low 
temperature the zero field resistivity,.dependence seems 
to be that of metallic ferromagnetstm. Also transport 
features at fixed magnetization m indicate that for large 
or small magnetization values the gxnonentially stretched 
dependence is markedly differentQErbil. The role of lat- 
tice strain as shown on films grown on different substrates 
does not seem to explain the findings, and should be sup- 
plemented by taking into account magnetic domain like 
effectsEi This implies that magnetic features should be 
reemphasized. The DES has been recently also much 
criticized because it was thought that it could not ex- 
plain the qualitative difference in electrical conduction 
for the whole range of x in one of the most often stud- 
ied CMR materials, i.e. La\^ x Ca x MnO^. However, it 
has been shown that a more correct treatment of the 
DES leads to an acceptable view in the interesting dop- 
ing range [x = 0.16,0.40], i.e. where the magnetic and 
conducting transitions are high (near room temperature). 
Outside this interval, further work should be envisaged 
since the magnetic structure has quite another period- 
icity, but the following concepts should still hold as it 
will be easily seen. In fact, three simple ingredients can 
be combined in a realistic way in order to emphasize the 
magnetic degrees of freedom role, i.e. the Ising model 
for the magnetic spins, the Drude formula for the electri- 
cal conductivity, and some type of scattering for hopping 
electrons. By combining these basic ingredients, we es- 
sentially take into account magnetic cluster effects, but 
the more so their correlated fluctuations. The former 
scattering strength and magnetic state life time are the 
only microscopic physical parameters. Notice that 

(i) Grains are usually pretty small, whence there is 
much grain boundary scattering. We neglect such a (ex- 
trinsic) background term here, though it can be easily 
inserted in the scattering strength if necessary. 

(ii) The number of carriers is also kept field and tem- 
perature independent. To take into account a density 
of states temperature (and even field) dependence is a 
rather trivial generalization to be made, within a self- 
consistent picture taking into account band and localized 
state carriers (with possible spin states). n 

(iii) The temperature dependent lattice distortioro role 
for hopping charge carriers is also neglected here, but 
can be included at a later stage again in the definition of 



1 



a temperature dependent lattice parameter and mapped 
into an effective carrier mass or into an effective localized 
spin coupling strength. 

In so doing, we do not claim that the present model 
should give precise quantitative values at this time be- 
cause of the extremely limited number of parameters 
which we are using. Nevertheless the theory will be 
in good agreement with experimental data. Moreover, 
extensions seem easy in light of the paraphernalia of 
solid state physics ideas and techniques. This model is 
surely not the unique alternative to the magnetic polaron 
model. However, this paper shows that the correlation 
between spins is the key ingredient to be integrated in 
the understanding of the CMR phenomenon, seen as a 
transport property in presence of magnetic states rather 
than a set of such near equilibrium specific states con- 
trolled by some unknown exchange interaction. For the 
following developments, let us reduce the problem to a 
two dimensional (d — 2) case, allowing as a first approx- 
imation an in-plane conduction like in thin filmscl. The 
Ising modelO on a d — 2 square lattice is used for the 
spins on the manganese sites assumed to represent the lo- 
cal magnetization of the system. This simplified picture 
allows for a faster way of obtaining the following results, 
but the spin-spin exchange interactions could be as well 
of indirect origin as in DES without loss of generality. 
Simply we let each lattice site i contains a two-state spin 
Ui = ±1. The dimcnsionless Hamiltonian reads 

E = -K^maj -h^ai (1) 

i,j i 

where K — J/kT is the dimensionless interaction be- 
tween nearest neighboring i,j spins, and h is the dimcn- 
sionless magnetic field H/kT orienting the spins. The 
(h = 0, K > 0) case is a classical problem taught in 
classrooms because it has a non-trivial phase transition 
as demonstrated by Onsager, i.e. a logarithmic diver- 
gence of the specific heat near the reduced critical tem- 
perature K c = — iln(\/2— 1). The Ising model implies 
that there are droplets (clusters) of, e.g. +1 or —1 spins 
which nucleate, grop. coalesce and disappear as a func- 
tion of temperatureEj. It is well knownE^rEj that a ferro- 
paramagnetic transition takes place exactly at K c for 
h = on such a lattice. In each grain or if the inter- 
grain coupling is adequate, clusters of respectively up 
and down spins coexist and the average size £ of these 
clustew diverges at K c following £ ~ \K — K c \^ lJ with 
v = lE-3. The other properties like m(K, h) are not known 
exactly because the Ising model in a field has not yet been 
solved. In view of the partially covalent-ionic bonding in 
the plane, the quasi localized carriers are supposed to 
be (spinless) electrons having a linear hopping motion 
along the electrical field imposed across the lattice. The 
Lorentz force is neglected here because of the rather short 
mean free path. In the computer experiments, we launch 
electrons toward the right at random from the left side 
of the "sample" . Each electron jumps to the right from a 



site to the next nearest neighbor site at each time step as 
follows. In a so-called "magnetic cluster", the electron 
hopping is free. When an electron reaches a magnetic 
cluster wall, the electron is stopped with a probability 
(1 — p) or transmitted with a probability p. By analogy 
with tunneling effect, p is assumed to be the exponen- 
tial of some measure of the cluster size s ahead of the 
electron, i.e. 

P = cxp (-7s) (2) 

where 7 is a dimensionless parameter which is like a po- 
tential barrier strength of the cluster. There is no reten- 
tion time upon a site nor phonon nor magnetic drag, nor 
other type of scattering. At each time step the magnetic 
structure is recalculated according to a Monte-Carlo pro- 
cedure for the Ising Hamiltonian. We count the carrier 
arrival time r on the right hand side of the "sample". 
This time obviously depends on the sign distribution fluc- 
tuations of the spins for a given K and h on the line 
during the electron hopping. 

Following the Drude formula, the resistivity is directly 
obtained from 

P=l (3) 

where L is the size of the lattice. At high temperature, 
when the spins are completely disordered the resistivity 
is of course large; it is smaller but not negligible at lower 
temperature; near the critical point K c for h — 0, the 
resistivity should become enormous: indeed the electron 
is a little bit "at a loss" because the spin fluctuations are 
huge and much hamper the electron motion. A magnetic 
field stiffens the clusters (or reduces the fluctuations). 
Therefore, the resistivity should be reduced because the 
electron has a greater chance to find its way through. 
Thus, the qualitative features of CMR are immediately 
found in this simple model. Notice that this CMR ver- 
sion is somewhat like a temperature dependent "mag- 
netic Kronig-Penney model" in an electric field since each 
wall is a potential barrier of which strength 7 is controlled 
by magnetic and therm aLconditions (Fig. lb), just like in 
disordered thin filmalj'EZI. The non-trivial (new) ingre- 
dient is that the "barriers" are correlated and controlled 
by the thermodynamics of the Ising model, in space and 
time. In order to obtain a good numerical convergence, 
we have left the magnetic system to reach a pseudo- 
steady-state before launching the "electrons". Lattices 
up to 256 x 256 were used. For conciseness, we fixed arbi- 
trary herein 7 = 1. The results do not change drastically 
with 7. Clearly, at a later stage of investigations, the 7s 
term can be itself temperature and/or magnetization de- 
pendent for taking into account the lattice strain. More 
complicated schemes taking into account different spin 
channels can be also imagined within an effective medium 
approximation. Here we consider that the spin density 
corresponds to the case where the majority of polarized 
(up or down) spins is much larger than the minority. The 
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spinless approximation of the charge carrier is not even 
a strong approximation. Indeed, it is clear that the true 
electron hopping is depending on the availability of a 
neighboring state of similar nature. Thus, the electron 
will be rather stopped in front of a wall after which the 
spins states have the opposite sign to that of the incom- 
ing electron. This is also in the DES spirit in fact. Figure 
2 presents the resistivity p as a function of K. Both zero 
and non-zero magnetic field cases are illustrated. The 
ferro-paramagnetic transition at K c is indicated by the 
vertical line. As expected, a bump is observed in p below 
K c and an inflexion at K c . Moreover, the bump height 
decreases as h increases. Such a bump can be viewed 
as the signature of a percolation transitionlia rather than 
a strict metal- insulator transition. On both sides of K c , 
the resistivity p decreases exponentially. This shows that 
the experimentally observed decreasing behavior in the 
high temperature phase might have nothing to do with 
a semiconducting phase or a "metal-insulator" transition 
as often claimed, but rather to the number of available 
final state in the scattering. The dimensionless excess 
resistivity defined by 



Ap = 



P(0) - P(h) 
P(0) 



(4) 



is shown as a function of K in Figure 3. Data due to 
different values of the magnetic field are shown. This 
quantity Ap is found to be independent of the length 
L of the lattice. The quantitative features of the CMR 
are well observed, i.e. (i) a peak at some intermediate 
temperature, (ii) an increase of Ap with h, and (hi) a 
shift of the peak towards high temperatures and (iv) a 
wide transition region. A non-trivial test of the model 
and theory is in order on Figure 4 which presents the 
evolution of the maximum of Ap as a function of ft. in a 
semi-log plot. The value of Ap for K = K c is also given. 
A logarithmic increase of Ap with h is observed, i.e. 



Ap ~ Inh. 



(5) 



Such a behavior can be found in La\- x Ca x Mru 
Lai- x Mg x MnOs, and Prx- x Sr x MnOz compounda§Ej 
(Figure 5). The predicted logarithmic behavior (straight 
lines) seems to hold quite well for moderate to high mag- 
netic fields. It is true that one should distinguish be- 
tween low field, moderate and high field CMREj. At 
lower field values, (e.g. H < 0.2 T) inhomogeneities play 
a relevant roleclcil, an effect outside the present investiga- 
tion. In addition, since CMR can be seen as related to a 
second-order critical phenomenon, it should be possible 
to describe CMR with the help of scaling arguments. For 
h = 0, the time for an electron to cross the system of size 
L is approximately given by 



L 



(6) 



i=l 



for £ < L where £ is the characteristic size of the clusters 
and (1 — exp (—7s)) 1 is the probability that the electron 



remains blocked i successive times on a magnetic wall be- 
fore being transmitted to the next site. Using the Drudc 
formula, one has 



1 



£ In (1- exp (- 7 0)' 



(7) 



The above relationship neglects the correlations between 
wall fluctuations. Nevertheless, the form of Eq.(8) ex- 
plains the finite size effects, (Figure 2), to be seen when 
£ w L and for small L values. Close to the ferro- 
paramagnetic transition, one can develop Eq.(8) and ob- 
tain 



p ~ exp (£) ~ exp 



1 



\K-K r 



(8) 



The wall fluctuation correlations being in fact negligeable 
far away from K c , the above exponential scaling behav- 
ior seems to be correct on both sides of K c . A quanti- 
tative comparison to available data is not immediate at 
low field (H < 300mT) since grain boundary background 
is very much field dependent there. In such a regime be- 
cause of inherent inhomogeneity of (usually small) grains, 
the Ap(hX dependence is indeed to be known as non 
universa£3. In the case of large h values, it has been 
argued elsewhereEj that the in-plane correlation length £ 
near K r scales as 



-u/0,5) 



(9) 



where 6=1/8 and S — 15 for the d = 2 Ising model 
whence the Ap ~ In h behavior should not be expected. 
Introducing the latter scaling relation in Eq.(9) leads to 
a stretched exponential, a law which qualitatively implies 
some apparent parallelism of Ap(h) decaying curves with 
an amplitude being temperature dependent. This feature 
can be understood as resulting from the non-negligable 
far away from K c cluster fluctuations and from some drift 
due to the magnetic field h. 

It may be recalled that an analogous treatment to ours 
occured-jn the pioneer theoretical work of Fisher and 
LangercJ using a d — 2 Ising model instead of a true d = 
3 model for describing the experimental results on the 
resistivity at ferromagnetic transitions. A divergence was 
predicted at T c while experimentally an inflexion (only) 
was seen. This paradoxical situation was found to be due 
to using a 2D rather than a 3D Ising model. Extensions 
towards a better agreement.to d = 3 cases should follow 
the same generalizing linesa in the future for our CMR 
model. 
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FIG. 1. Schematic illustration of the magnetotransport 
process as discussed in the text: (a) The square lattice on 
which +1 and -1 magnetic domains are distributed. On this 
lattice, one electron follows the linear motion illustrated by 
the dashed arrow; (b) The barrier landscape viewed by this 
electron at the time corresponding to the snapshot of (a). 
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FIG. 2. Theoretical resistivity p as a function of K . Two plG 3 Magnetoresistance Ap as a f unc ti on of K. 

cases are shown: h — and h ^ 0. Different lattice sizes are 
illustrated: L = 32, L = 64. 
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FIG. 4. Semi-log plot of the evolution of the maximum 
of the magnetoresistance Ap as well as Ap for K = K c as a 
function of the magnetic field h. 



FIG. 5. Semi-log plot of the evolution of the magne- 
toresistance Ap as a function of the magnetic field for var- 
ious .compounds: Lai- x Ca x MnQa. fromU, La\- X M g x MnO-i 
fromfl and Pri- x Sr x MnC>3 fromEj. The logarithmic behav- 
ior (straight line) predicted by the model (Eq.(6)) is clearly 
observed. 



